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ABSTRACT 

We study the statistics of large-separation gravitational lens systems produced by non-spherical halos 
in the Cold Dark Matter (CDM) model. Specifically, we examine how the triaxiality of CDM halos 
affects the overall lensing probabilities and the relative numbers of different image configurations 
(double, quadruple, and naked cusp lenses). We find that triaxiality significantly enhances lensing 
probabilities by a factor of ~2— 4, so it cannot be ignored. If CDM halos have central density slopes 
a < 1.5, we predict that a significant fraction (>20%) of large-separation lenses should have naked 
cusp image configurations; this contrasts with lensing by isothermal (a ~ 2) galaxies where naked 
cusp configurations are rare. The image multiplicities depend strongly on the inner density slope a: 
for a = 1, the naked cusp fraction is >60%; while for a = 1.5, quadruple lenses are actually the most 
probable. Thus, the image multiplicities in large-separation lenses offer a simple new probe of the 
inner density profiles of dark matter halos. We also compute the expected probabilities and image 
multiplicities for lensed quasars in the Sloan Digital Sky Survey, and argue that the recent discovery 
of the large-separation quadruple lens SDSS J1004+4112 is consistent with expectations for CDM. 

Subject headings: cosmology: theory — dark matter — galaxies: clusters: general — gravitational 
lensing 



1. INTRODUCTION 

The Cold Dark Matter (CDM) model of structure for- 
mation naturally predicts the existence of strong gravi- 
tational lens systems with image separations of ~ 10" or 
even larger. Observations of massive clusters of galaxies 
have revealed many systems of "giant arcs" representing 
lensed images of background galaxies (Lynds & Petrosian 
1986; Soucail et al. 1987; Luppino et al. 1999; Gladders 
et al. 2003; Zaritsky & Gonzalez 2003). However, until 
recently all lensed quasars and radio sources had image 
separations < 7" corresponding to lensing by galaxies, 
despite some explicit searches for lenses with larger sep- 
arations (Phillips et al. 2001; Ofek et al. 2001). 4 Lensing 
of quasars by clusters was finally observed with the re- 
cent discovery and confirmation of SDSS J1004+4112, a 
quadruple lens with an image separation of 14'.'62 found 
in the Sloan Digital Sky Survey (Inada et al. 2003; Oguri 
et al. 2004). This lens confirms an important predic- 
tion of the CDM model; indeed, the lensing probability 
inferred from the discovery is in agreement with reason- 
able values of the cosmological parameters (Oguri et al. 
2004). 

The statistics of large-separation lenses can be used 
to place constraints on the density profile of dark ha- 
los (Maoz et al. 1997; Keeton & Madau 2001; Keeton 
2001a; Wyithe, Turner, & Spergel 2001; Takahashi & 
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4 Miller et al. (2004) recently reported six pairs of quasars in 
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but none of the candidates has been confirmed, and theoretical ar- 
guments by Oguri (2003) indicate that it would be quite surprising 
if any of the systems are lenses. 



Chiba 2001; Sarbu, Rusin, & Ma 2001; Li & Ostriker 
2002; Oguri et al. 2002; Oguri 2002; Huterer & Ma 2004; 
Kuhlen, Keeton, & Madau 2004), or determine the abun- 
dance of massive dark halos (Narayan & White 1988; 
Wambsganss et al. 1995; Kochanek 1995; Nakamura & 
Suto 1997; Mortlock & Webster 2000; Oguri 2003; Lopes 
& Miller 2004; Chen 2004). The discovery of the first 
such lens suggests that these statistics can be a practical 
tool to study structure formation in the universe. The 
statistics of giant arcs are also known as a good probe of 
clusters (Bartelmann et al. 1998; Meneghetti et al. 2001; 
Molikawa & Hattori 2001; Oguri, Taruya, & Suto 2001; 
Oguri, Lee, & Suto 2003; Wambsganss, Bode, & Ostriker 
2004; Dalai, Holder, & Hennawi 2004; Maccio 2004), and 
in fact lensed arcs and quasars complement each other in 
several ways. For instance, in lensed quasar surveys one 
first identifies source quasars and then checks whether 
they are lensed, while in searching for lensed arcs one 
selects massive clusters and then searches for lensed arcs 
in them. In other words, surveys for arcs are biased 
toward high mass concentrations, while lensed quasars 
probe random lines of sight. Clusters selected by the 
presence of lensed quasars could, in principle, differ from 
those selected as having giant arcs. In addition, lensed 
quasars have three advantages over lensed arcs in sta- 
tistical studies. First, quasars can be regarded as point 
sources, while sources for arcs are galaxies whose intrinsic 
sizes and shapes are important but unobservable. Sec- 
ond, the number and configuration of images in a quasar 
lens system is unambiguous. Third, the redshift distri- 
bution of arc sources is poorly known (and controversial; 
see Oguri et al. 2003; Wambsganss et al. 2004; Dalai et 
al. 2004), while the redshift distribution of quasars is well 
known. 

In all previous analytic work on the statistics of large- 
separation lensed quasars, the lens objects were assumed 
to be spherical. However, in the CDM model dark ha- 
los are not spherical at all but triaxial (e.g, Jing & Suto 
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2002, hereafter JS02). It is already known that triaxial- 
ity has a significant effect on the statistics of lensed arcs, 
from both analytic (Oguri ct al. 2003, hereafter OLS03) 
and numerical (Meneghetti, Bartelmann, & Moscardini 
2003a; Dalai et al. 2004) points of view. In the statistics 
of normal lensed quasars, triaxiality (or ellipticity) has 
been thought to mainly affect the image multiplicities, 
with only small changes to the total lensing probabil- 
ity (Kochanek 1996; Keeton, Kochanek, & Seljak 1997; 
Evans & Hunter 2002; Chae 2003; Huterer & Keeton 
2004). However, that conclusion is based on nearly- 
singular isothermal lens models, and the situation may 
be quite different for the less concentrated mass distri- 
butions of the massive halos that create large-separation 
lenses. Moreover, only triaxial modeling allows us to 
study image multiplicities, and to consider whether it is 
statistically natural that the first known large-separation 
lens is a quadruple. 

The structure of this paper is as follows. In §2 we 
review the triaxial dark halo model and its lensing prop- 
erties (originally presented by JS02 and OLS03). In §3 
we show our general results, while in §4 we customize 
our predictions to the SDSS quasar sample. We sum- 
marize our conclusions in §5. Throughout the paper, we 
assume a A-dominated cosmology with current matter 
density CIm = 0.3, cosmological constant = 0.7, di- 
mensionless Hubble constant h = 0.7, and normalization 
of matter density fluctuations a$ — 0.9. 

2. FORMALISM 
2.1. Lensing by triaxial dark halos 

In this section, we briefly summarize the lensing prop- 
erties of the triaxial model of dark halos proposed by 
JS02. For more details, please refer to JS02 and OLS03. 

First, we relate the principal coordinate system of the 
triaxial dark halo x — (x, y, z) to the observer's coordi- 
nate system x' — (x', y' , z'), where the z'-axis runs along 
the line of sight to the observer. In general, the coordi- 
nate transformation is expressed as x — Ax' with 
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The density profiles of triaxial dark matter halos pro- 
posed by JS02 (also see Zhao 1996) is 



p(R) = 
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where 
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Two models commonly discussed in the context of CDM 
simulations are a — 1 and 1.5 (Navarro, Frenk, & White 
1997; Fukushige & Makino 1997; Moore et al. 1999; Jing 
& Suto 2000; Power et al. 2003; Fukushige, Kawai, & 
Makino 2004), and we focus on these. JS02 give fitting 
formulas for the axis ratios a/c and a/b in the triaxial 
model, and for the concentration parameter c e = R e / Ro, 
where R e is defined such that the mean density within 
the ellipsoid of the major axis radius R e is A e Q(z)p a i t (z) 

with A e = 5A vir (c 2 /ab)°- 75 . 



What matters for lensing is the projected surface mass 
density in units of the critical density for lensing, or the 
convergence k, which can be expressed as (OLS03) 
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where &tnfw is a dimensionless "strength" parameter 
(defined in OLS03), q = q y /q x < 1 is the axis ratio of the 
projected mass distribution, and 
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For a — 1, equation (5) has an analytic expression 
(Bartelmann 1996). For a = 1.5, we adopt a fitting for- 
mula for equation (5) presented by OLS03. Note that the 
convergence has elliptical symmetry, so the lensing de- 
flection and magnification can be computed with a set of 
1-dimensional integrals (Schramm 1990; Keeton 2001b). 
Also note that if we work in dimensionless units, scaling 
all lengths by Lq = Roq x , then the lensing properties of 
the dark matter halos depend only on the parameters a, 
^TNFWi and q. 

2.2. Cross sections and image separation distributions 

We compute lensing cross sections using Monte Carlo 
methods. Working in dimensionless coordinates X = 
x' I Lq and Y = y' /Lq, we pick random sources and use 
the gravlens software by Keeton (2001b) to solve the lens 
equation. Figure 1 shows examples of the three different 
kinds of image configurations: double, quadruple, and 
naked cusp lenses. 5 We count the number of sources that 
produce lenses of different image multiplicities to deter- 
mine the dimensionless cross sections 02, 04, and a c for 
doubles, quadruples, and cusps, respectively. For each 
set of images, we define the dimensionless image separa- 
tion 9 to be the maximum separation between any pair of 
images; this is a convenient definition that depends only 
on observable quantities and is well defined for all image 
configurations (no matter how many images there are). 
We bin the sources by the image separations they pro- 
duce to derive image separation distributions, as shown 
in Figures 2 and 3. For a given halo there is a range of 
separations, but it tends to be fairly narrow (<20%); the 
main exception is for cusp configurations, which show a 
tail to small separations that corresponds to sources near 
the cusp in the caustic. 

If we not only count the sources but also weight them 
appropriately, we can compute the magnification bias. 
Specifically, if the sources have a simple power law lu- 
minosity function 4>l(L) cx L~P then the "biased cross 
section" can be written as 



-I 



dXdY 



cj> L {L/p)/p 
4>l{L) 



= j dXdY^- 1 , 



(6) 



where the integral is over the multiply-imaged region of 
the source plane. We can compute the biased cross sec- 
tions for doubles, quadruples, and cusps similarly. Each 

5 We use the terms "double" and "quadruple" because the third 
and fifth images arc usually too faint to be observed, although with 
the density profiles we use here they are probably not as faint as 
for nearly-isothermal lenses (see Rusin 2002). 
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Fig. 1. — Sample image configurations. The top panels show the source planes, and the bottom panels show the corresponding image 
planes. The solid lines indicate the caustics and critical curves. We show three sources (denoted by triangles, circles, crosses), and their 
corresponding images. From left to right, the lenses are doubles, quadruples, and cusps. Specific values of (&tnfw s O) f° r each example 
are (2,0.95), (2,0.75), and (0.6,0.25) for doubles, quadruples, and cusps, respectively. Doubles and cusps are distinguished by the image 
parities: doubles have one positive-parity image and one negative-parity image, plus a central double-negative image that is usually too 
faint to be observed; while cusps have two positive-parity images and one negative parity image, all of comparable brightnesses. 



source is to be weighted by n 13 ^ 1 , where we take [i to 
be the magnification of the second brightest image to re- 
flect the method of searching for large-separation lenses 
in observational data such as the SDSS (see Inada et al. 
2003; Oguri et al. 2004). 

An important qualitative result is already apparent 
from Figure 1. CDM-type dark matter halos are very 
sensitive to departures from spherical symmetry, in the 
sense that even small projected ellipticities lead to large 
tangential caustics and hence large quadruple cross sec- 
tions. When the ellipticity is large, the tangential caustic 
is much larger than the radial caustic and nearly all the 
images correspond to cusp configurations. This situa- 
tion is notably different from what happens in lensing by 
galaxies that have concentrated, roughly isothermal mass 
distributions. In that case, the ellipticity must approach 
unity before cusp configurations become common (see 
Keeton et al. 1997; Rusin & Tegmark 2001). Such large 
ellipticities are uncommon, and cusp configurations are 
correspondingly rare among observed galaxy-scale lenses: 
among ^80 known lenses there is only one candidate 
(APM 08279+5255; Lewis et al. 2002). The incidence of 
cusp configurations therefore appears to be a significant 
distinction between normal and large-separation lenses. 



2.3. Lensing probabilities 

The probability that a source at redshift zs is lensed 
into a system with image separation is computed by 
summing the biased cross section over an appropriate 
population of lens halos: 



dP 
~d6 



x J d(a/b) J old J d(j) 



(1 + z L ) 3 / d(a/c) J dc e 
p(a/c)p(c e ) 



rffl 

p(a/b\a/c)p(e)p(4,)Ba — 



•(7) 



M(6) 



The first integral is over the volume between the ob- 
server and the source. The next three integrals are over 
the structural parameters of the lens halos, while the last 
two integrals cover the different orientations. The mass 
function of dark matter halos is dn/dM, and we use the 
model from equation (B3) of Jenkins et al. (2001). Fi- 
nally, M{9) is the mass of a halo that produces image 
separation 9 (for given redshift and other parameters), 
which is given by the solution of 



6 — R q x $(&TNFWj 



(8) 



The square brackets in equation (7) indicate that the 
integrand is to be evaluated only for parameter sets that 
produce the desired image separation. Equation (7) gives 
the total lensing probability, but we can simply replace 
the total biased cross section Bo with B 2 o- 2 , B404, or 
B c o c to compute the probability for doubles, quadruples, 
or cusps. 

We use the model given by JS02 for the probability 
distribution functions (PDFs) that appear in equation 
(7). The PDFs for the axis ratios are: 
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Fig. 2. — Image separation distributions for sample lenses with a = 1. Arrows indicate the average separations. The corresponding 
caustics are shown for reference. For each 6tnfWi the caustics are all plotted on the same scale. 




^ /2a/b — 1 — max(a/c, 0.5) 
\ 1 — max(a/c, 0.5) 



,(10) 



where M* is the characteristic nonlinear mass such that 



the RMS top-hat-smoothed overdensity at that mass 
scale is 1.68. Note that M* and fl(z) depend on red- 
shift, so p(a/c) varies with redshift such that halos tend 
to be more triaxial at higher redshifts (see Fig. 8 of JS02). 
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Fig. 4. — Lensing probabilities and image multiplicities as a 
function of a/c. The left panels are for inner slope a = 1, and the 
right panels for a = 1.5. We adopt 8 = 15", (i = 2.5, and 2:3 = 2.0. 
We show the total lensing probability (thick solid line), as well 
as the probabilities for double (dotted), quadruple (dashed), and 
cusp (dash-dotted) lenses. For comparison, the lensing probability 
of spherical halos is shown by the thin solid line. 



Also note that p(a/b\a/c) = for a/b < max(a/c, 0.5). 
The PDF for the concentration parameter is 
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For the median concentration parameter c e , we adopt a 
fitting formula given by JS02 (see also OLS03): 



c e = 1.35 cxp 



0.3 



'Avir(z c ) fl 
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with z c being the collapse redshift of the halo of mass 
M v ; r and A e = 1.1 in the Lambda-dominated CDM 
model (JS02). For spherical halos, large-separation lens 
statistics are highly sensitive to the concentration dis- 
tribution (both the median and the scatter; see Kuhlcn 
et al. 2004). However, we checked that in the triaxial 
model the sensitivity is much reduced because there is 
such a broad distribution of axis ratios; hence we do not 
present results with a different median or scatter in the 
concentration distribution. Finally, the PDFs for the ori- 
entation angles are 

sin6> 



1 

2V 



(13) 
(14) 



corresponding to random 3-d orientations. 



3. LENSING PROBABILITIES AND IMAGE 
MULTIPLICITIES IN THE TRIAXIAL HALO MODEL 

3.1. Dependence of the triaxiality 

We begin by examining how the lensing probabilities 
and image multiplicities vary when we change the de- 
gree of triaxiality. We remove the integral over a/c in 
equation (7) to compute the lensing probabilities at fixed 




£ 10 



100 10 
0[arcsec] 

Fig. 5. — Lensing probabilities and image multiplicities with 
triaxial dark halos as a function of image separation 9. The source 
is placed at 23 = 2.0, and slope of the source luminosity function 
is fixed to /3 = 2.5. 



triaxiality (we still integrate over the intermediate axis 
ratio a/b and over random orientations). We can then 
plot the probabilities as a function of a/c, as shown in 
Figure 4. In this example, we place the source at red- 
shift z$ = 2.0, and we use a source luminosity function 
with slope (3 = 2.5. We compute the probabilities for an 
image separation of 8 = 15", similar to that of the one 
known large-separation lensed quasar SDSS J1004+4112 
(Inada et al. 2003). 

For a/c -^lwc recover the spherical case. As a/c de- 
creases (the triaxiality increases), at first the total lens- 
ing probability stays roughly constant but the fraction of 
quadruples rises; this is similar to the effects of elliptic- 
ity on isothermal lenses (see Keeton et al. 1997; Rusin & 
Tegmark 2001). Then the probability for naked cusp im- 
age configurations begins to rise dramatically, and they 
come to dominate the total probability. Interestingly, the 
sum of the probabilities for quadruple and double lenses 
is roughly equal to the probability for spherical halos for 
most values of a/c, at least for this example with image 
separation 9 = 15" (also see Figures 5-7 below). This 
suggests that the enhancement in the total lensing proba- 
bility is mainly driven by naked cusp configurations. The 
a = 1 and 1.5 cases both have these qualitative features, 
and they differ only in the quantitative details. Since 
the typical triaxiality in CDM simulations is a/c ~ 0.5 
(JS02), it appears that triaxiality can have a significant 
effect on the statistics of large-separation lenses. 

3.2. Full results 

To compute the full impact of triaxiality on lens statis- 
tics, we must integrate over an appropriate triaxiality 
distribution (as in equation 7). Figure 5 shows the re- 
sulting lensing probabilities and image multiplicities as 
a function of the image separation 9. Again, we place 
the source quasar at zs = 2.0, and fix the slope of the 
source luminosity function to (3 = 2.5. The first im- 
portant result is that the triaxial model predicts larger 
lensing probabilities than the spherical model for all im- 
age separations. The enhancement is a factor of ~4 for 
a = 1, and a factor of ^2 for a = 1.5, if the image sep- 
aration is not so large (9 < 30"). At larger separations 
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Fig. 6. — Lensing probabilities and image multiplicities with tri- 
axial dark halos as a function of the slope of the source luminosity 
function (3. We consider an image separation of 8 = 15", and we 
place the source at zg = 2.0. 
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Fig. 7. — Lensing probabilities and image multiplicities as a 
function of source redshift zg. We fix the image separation to 
8 = 15" and the of the source luminosity function to fi = 2.5. 



it seems that the effect of triaxiality is even more signif- 
icant, especially for a = 1; we will discuss this issue in 
§3-3. 

There are several interesting results in the image mul- 
tiplicities. The a = 1 and 1.5 cases have very different 
multiplicities: with a = 1 the lensing probability is domi- 
nated by cusp configurations; while with a — 1.5 quadru- 
ple lenses are somewhat more common than cusps. Nei- 
ther result is very sensitive to the image separation. In 
both cases double lenses are fairly uncommon, which is 
very different from the situation with normal arcsecond- 
scale lenses produced by nearly-isothermal galaxies. This 
result is consistent with previous theoretical conclusions 
that image multiplicities depend on the central concen- 
tration of the lens mass distribution, such that less con- 
centrated profiles tend to produce more quadruple and 
cusp lenses (Kassiola & Kovner 1993; Kormann, Schnei- 
der, & Bartelmann 1994; Rusin & Tegmark 2001; Evans 
& Hunter 2002: Dalai et al. 2004). The important point 
for observations is that if dark halos have inner profiles 
with a < 1.5, then many or even most large-separation 
lenses should be quadruples or cusps rather than doubles. 
(We will consider the implications for SDSS J1004+4112 
in §4.) Another point is that the image multiplicities are 
sensitive to the inner density profile, so they offer a new 
method for probing dark matter density profiles that is 
qualitatively different from methods discussed before. 

Magnification bias is important in lens statistics, par- 
ticularly in image multiplicities, because it gives more 
weight to quadruple and cusp configurations (which tend 
to have large magnifications) than to doubles. We should 
therefore understand what happens when we modify the 
magnification bias by varying the slope (3 of the source lu- 
minosity function. The results are shown in Figure 6 (for 
image separation 9 = 15" and source redshift z$ = 2.0). 
The lensing probabilities increase as (3 increases, because 
as the source luminosity function becomes steeper mag- 
nification bias becomes stronger. 6 Interestingly, the in- 
crease in the total probability due to triaxiality weakens 

6 Note that the magnification bias diverges if the luminosity 
function is a pure power law with fi > 3, so we are restricted to 
shallower cases. 



as (3 increases, although the effect is not strong. As for 
the image multiplicities, a = 1 halos are always domi- 
nated by cusp lenses, although for sufficiently steep lu- 
minosity functions quadruples become fairly common. 
With a = 1.5 halos, when magnification bias is weak 
(/3 <~ 1) doubles are the most probable, but as magnifi- 
cation bias strengthens (/3 increases) quadruples receive 
more weight and become the most likely. In practice, the 
effective values of j3 are larger than ~1.5 for both optical 
(e.g., Boyle ct al. 2000) and radio surveys (e.g., Rusin & 
Tegmark 2001), so we expect cusps to dominate for a = 1 
and quadruples to be the most common for a = 1.5. 

Finally, we consider whether the results depend on 
the source redshift, as shown in Figure 7. The lensing 
probabilities rise with zs, because there is more volume 
and hence more deflectors between the observer and the 
source. In addition, geometric effects mean that a given 
halo can produce a larger image separation when the 
source is more distant, so the halo mass required to pro- 
duce a given image separation goes down and the abun- 
dance of relevant deflectors goes up. However, the prob- 
ability increase affects the different image configurations 
in basically the same way, so the image multiplicities are 
quite insensitive to the source redshift. Therefore, we 
conclude that details of the source redshift distribution 
are not so important for image multiplicities, at least 
when the source luminosity function has a power law 
shape. 

3.3. Statistics at larger image separations 

In Figure 5, at very large image separations (9 > 100") 
the lensing probabilities in the triaxial halo model are or- 
ders of magnitude larger than those in the spherical halo 
model. In addition, at these very large separations the 
a = 1 case produces higher probabilities than the more 
concentrated a = 1.5 case. Both features are puzzling 
and invite careful consideration. 

Figure 8 shows the dependence of total lensing proba- 
bility on the lower limit of the integral over a/c. In the 
previous calculations, we assumed (a/c) m ; n = 0.1. This 
figure shows that for 9 — 15" the results are quite insen- 
sitive to (a/c) m i n , suggesting that the contribution from 
extremely triaxial halos is negligible. For 9 — 200", how- 
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Fig. 8. — Dependence of the lensing probability on the cutoff 
in a/c. Both 9 = 15" and 200" are shown. Filled squares and 
open circles denote a = 1 and 1.5, respectively. The slope of the 
luminosity function is fixed to (3 = 2.5. 



ever, the lensing probability rapidly decreases as (a/c) m i n 
increases. In other words, the lensing probability at very 
large image separations seems to be dominated by very 
small a/c, or very large triaxialities. 

Results that are dominated by such extreme halos are 
probably not very reliable. They depend sensitively on 
both the assumed PDF for the axis ratio a/c (eq. [9]) 
and the correlation between a/c and the concentration 
c e (cq. [12]) at very small axis ratios. The fitting forms 
presented by JS02 were intended to reproduce the PDF 
and correlation at a/c > 0.3, and it is unclear whether 
they are still accurate at a/c ~ 0.1. In addition, even if 
we know accurate fitting forms, such a situation implies 
that sample variance (i.e., the effect of the finite number 
of lensing clusters) may be quite large. 

Another ambiguity is the projection effect. In this pa- 
per, we assumed that the density profile (eq. [2]) extends 
beyond the virial radius, and in projecting along the line 
of sight we integrated the profile to infinity. Although 
it is not clear whether we should cut off the profile at 
the virial radius or not (e.g., Takada & Jain 2003), the 
effect of the extended profile on the gravitational lensing 
is not so large for normal dark halos. However, when a/c 
is small enough, equation (12) indicates that the con- 
centration parameter c e becomes smaller than unity, so 
the effect of the extended profile outside the virial ra- 
dius is quite significant. The puzzling feature that the 
a = 1 case produces higher probabilities than the a = 1.5 
case at very large separations can be ascribed to the pro- 
jection effect, because the effect is more significant for 
shallower density profiles. 

Thus, lenses with extremely large image separations 
are associated with the most extreme dark matter halos, 
and it may be difficult to make reliable predictions about 
them. We emphasize, though, that these issues do not 
apply to lenses with separations 9 < 30", and on these 
scales we believe our results to be robust. 

4. PREDICTIONS FOR THE SDSS 

In the previous section, we assumed a simple power law 
source luminosity function to understand the general ef- 
fects of triaxiality. Here we consider a specific luminosity 



function appropriate for quasars, and further customize 
our predictions to quasars in the Sloan Digital Sky Sur- 
vey. 

We adopt the double power law 5-band luminosity 
function for quasars proposed by Boyle, Shanks, & Pe- 
terson (1988), 



b L (z s ,L)dL = 



dL 



[L/L.{z s )]* + [L/L.{z 8 )]& L.(zs)' 

(15) 

and use (3 h = 3.43 and fa = 1.64 (see Boyle et al. 2000). 
We let the break luminosity evolve following a model that 
reproduces the low-redshift luminosity function as well as 
the space density of high-redshift quasars (see Wyithe & 
Loeb 2002; Oguri et al. 2004). In practice we actually 
use the cumulative luminosity function 



-i: 



b L (z s ,L)dL, 



to calculate the biased cross section (see §2.2) 



Bb 



-I 



dXdY- 



(16) 



(17) 



We approximate that the SDSS quasar sample is a sam- 
ple with a flux limit of i* = 19. 7. 7 One needs a cross-filter 
K-correction to convert observed i* magnitudes to abso- 
lute i3-band luminosity. We adopt the following form 



K Bi (z) = -2.5(1 -a s )log(l + z) 



- 2 ' 5asl0g 4400 



- 0.12, 



(18) 



with a s = 0.5 (Oguri et al. 2004). Finally, we approxi- 
mate the redshift distribution (sec Fig. 1 of Oguri et al. 
2004) with the following Gaussian distribution: 

: - 1.45) 2 



p(z s )dz s 



1 

L21 



exp 



dz s 



2(0.55) 2 

(0.6 <z s < 2.3). (19) 

We have confirmed that the results using these approx- 
imations agree well with those obtained by fully taking 
account of the observed redshift and magnitude distribu- 
tions. 

Figure 9 shows the lensing probabilities and image mul- 
tiplicities as a function of image separation 9 for the 
SDSS quasar sample. Again large-separation lenses are 
dominated by cusp configurations for a = 1, but all 
three configurations are almost equally likely for a = 1.5. 
Therefore, we confirm that image multiplicities in SDSS 
large-separation lenses will offer interesting information 
on the density profile of dark halos. We can now consider 
whether it is statistically natural that the first large- 
separation lens in the SDSS is a quadruple lens. We 
find that for an image separation 9 = 15" the fractions 
of quadruple lenses are ~0.2 and ~0.4 for a = 1 and 1.5, 
respectively. Thus a = 1.5 could explain the discovery 
of the quadruple lens somewhat better, but a = 1 is also 
not unnatural. 

7 The SDSS quasar target selection is aimed to choose quasars 
with i* < 19.1 (Richards et al. 2002). However, we assume a flux 
limit of i* = 19.7 because there are quasars with i* > 19.1 in 
the SDSS quasar sample which are, for example, first targeted as 
different objects but revealed to be quasars. 
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Fig. 9. — Lensing probabilities and image multiplicities for SDSS 
quasars at rcdshifts 0.6 < z$ < 2.3. 



Finally, we can use the discovery of SDSS J1004+4112, 
together with the lack of large-separation lenses in the 
Cosmic Lens All-Sky Survey (CLASS; Phillips et al. 
2001), to constrain the cosmological parameter as de- 
scribing the normalization of the density fluctuation 
power spectrum. Using spherical models, Oguri et al. 
(2004) found that the discovery of SDSS J1004+4112 re- 
quired rather large values of either a or as, but given 
the importance of triaxiality we should revisit this ques- 
tion. For the SDSS, we compute the expected number 
of large-separation lenses with 7" < 8 < 60" among the 
29811 SDSS quasars. For CLASS, we adopt a power law 
source luminosity function with (3 = 2.1 (see Rusin & 
Tegmark 2001), fix the source redshift to z$ = 1.3 (see 
Marlow et al. 2000), and calculate the expected number 
of lenses with 6" < 9 < 15" among 9284 flat-spectrum 
radio sources (Phillips et al. 2001). We then compute 
the likelihood 

C oc (1 - e - NsDSS ) e ~ Nc ^ ASS , (20) 

which represents the Poisson probability of observing no 
large-separation lenses in CLASS when Aclass are ex- 
pected, and at least one large-separation lens in SDSS 
when Asdss are expected. (There may be other large- 
separation lenses in the SDSS sample that have not yet 
been identified.) We consider two possibilities for the ex- 
pected number of lenses in the SDSS: (1) the total num- 
ber of lenses A tot is used as Asdss; (2) only the number 
of quadruple lenses A qua d is because the discovered lens 
is quadruple. 

Figure 10 shows the resulting maximum likelihood con- 
straints on as- We find that as ~ 1 explains the data 
well, although the details depend on the value of a and 
the choice of Asdss- Only the case with a = 1 and 
Asdss = Aq Uad prefers relatively large a$, but a s = 1 
is still allowed at the 2a level. At present the data do 
not allow particularly strong constraints on as- Never- 
theless, we can conclude that the status large-separation 
lenses is quite consistent with the predictions of CDM 
given as ~ 1. 

For comparison, Figure 10 also shows results for spher- 
ical halos. It turns out that the spherical model overes- 
timates the value of a s by ^0.1 for a = 1 and ^0.2 for 
a = 1.5, compared with cases where we take Asdss = 



1 — i — i — i — i — | — i — i — i — i — |- 




FlG. 10. — Maximum likelihood estimates for ag from combin- 
ing the discovery of SDSS J1004+4112 in SDSS with the lack of 
large-separation lenses in CLASS. In making predictions for SDSS, 
we consider two cases: the appropriate prediction could be the to- 
tal number of lenses (solid lines); or since SDSS J1004+4112 is a 
quad the appropriate quantity could be the number of quadruples 
(dashed lines). The likelihoods for a = 1 and 1.5 are shown by 
thick and thin lines, respectively. Results for spherical halos are 
also shown by dash-dotted lines for reference. 



At t • Triaxiality is therefore an important systematic ef- 
fect in these cases. Interestingly, the best-fit values of 
as from spherical models are quite similar to those from 
triaxial models with Asdss = A quac j. In both cases, the 
likelihood function for the spherical model is narrower 
than for the triaxial model, indicating that the spherical 
model would underestimate the statistical uncertainties 
in a 8 . 

5. SUMMARY AND DISCUSSION 

The dark matter halos predicted by the CDM model 
are triaxial rather than spherical, which has a significant 
effect on the statistics of large-separation gravitational 
lenses. Triaxiality systematically enhances the lensing 
probability by a factor of ^4 if dark halos have an inner 
density profile with a = 1, or a factor of ^2 if a = 1.5. 
The effects may be even more dramatic at very large 
image separations {0 > 100"), although such lenses are 
very sensitive to the most triaxial halos and so the pre- 
dictions are not as reliable. Thus, triaxiality must be 
added to the list of important systematic effects that 
need to be included in calculations of large-separation 
lens statistics. (Some of the other effects are the inner 
density profile and the shape of the distribution of con- 
centration parameters, as found by previous studies of 
large-separation lens statistics using spherical halos.) 

Triaxial modeling allows us to predict the image mul- 
tiplicities for large-separation lenses. We found that the 
multiplicities depend strongly on the density profile: for 
a = 1, lenses are dominated by naked cusp image con- 
figurations; while for a — 1.5, quadruple configurations 
are the most probable. Double lenses, which are domi- 
nant among normal arcsecond-scale lenses, are subdom- 
inant in both cases. Note that cusp lenses can be distin- 
guished from doubles by the presence of a third image 
comparable in brightness to the other two, and by the 
configuration of image positions. The differences can be 
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ascribed to the different mass density profiles, and they 
indicate that the multiplicities of large-separation lenses 
will provide a qualitatively new probe of the central den- 
sity profiles of massive dark matter halos (and hence a 
new test of CDM). 

We have computed lensing probabilities and image 
multiplicities for the SDSS quasar sample. We predict 
that for both a = 1 and 1.5 most of the large-separation 
lenses should be quadruples or cusps. The fractions of 
quadruple lenses at separations of 6 = 15" are ^0.2 and 
~0.4 for a = 1 and 1.5, respectively. Thus it is not sur- 
prising that the first large-separation lens discovered is a 
quadruple. In addition, we computed the expected num- 
ber of large-separation lenses in both SDSS and CLASS, 
and found that the data are consistent with the CDM 
model with <r 8 ~ 1, in agreement with other measure- 
ments (e.g., Spergel et al. 2003, and references therein). 
Thus, the discovery of SDSS J1004+4112 can be inter- 
preted as additional support for CDM on non-linear clus- 
ter scales. 

The prediction that triaxial halos produce significant 
fractions of quadruple and cusp lenses should be kept 
in mind when considering samples of candidate large- 
separation lenses. For example, Miller et al. (2004) found 
six large-separation double lens candidates in the 2dF 
Quasar Redshift Survey, but no quadruple or cusp lens 
candidates. Even accounting for small number statistics, 
our results suggest that such a high fraction of doubles 
would be inconsistent with CDM at more than 3<7, 8 and 
that it would be surprising if many of the six candidates 
are genuine lens systems. 

We note that the triaxial dark halo model we adopted 
in this paper can be improved in several ways. First, we 
assumed that the axis ratios of the triaxial ellipsoids are 
constant with radius. However, JS02 showed that the 
axis ratios decrease slightly toward the halo centers: a/c 
decreases by ^0.2 as the mean radius decreases from 
^0.6r v ; r to ^0.06r v j r . Since strong lensing is most sen- 
sitive to the inner parts of dark halos, it is possible that 
we have actually underestimated the effects of triaxial- 
ity on the statistics of large-separation lenses. On the 
other hand, including baryons (which were neglected in 
the simulations of JS02) would tend to make dark halos 
rounder by A(a/c) <~ 0.1-0.2 because of the isotropic gas 
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pressure. 

While our theoretical model is much more realistic 
than the simple spherical model, we have still made sev- 
eral simplifying assumptions. One is that we have ne- 
glected substructure in dark halos. The galaxies in mas- 
sive cluster halos do not have a large effect on the statis- 
tics of lensed arcs (Meneghetti et al. 2000), but it is not 
obvious whether or not they would affect large-separation 
lenses. Substructure can affect the image multiplicities 
for isothermal lenses (Cohn & Kochanek 2004), so it 
should be considered for CDM halos as well. Another ef- 
fect we have neglected is the presence of a massive central 
galaxy in a cluster. Meneghetti, Bartelmann, & Moscar- 
dini (2003b) claim that central galaxies do not have a 
large effect on arc statistics. However, because our re- 
sults depend on the inner slope of the density profile, and 
a central galaxy effectively increases the concentration, 
this effect should be considered. A third phenomenon we 
have neglected is cluster merger events. Indeed, mergers 
can change the shapes of critical curves and caustics sub- 
stantially, and thus have a great impact on lensing cross 
sections (Torri et al. 2004). To estimate the effect on 
large-separation lens statistics, we would need a realistic 
model of the cluster merger event rate and the physical 
conditions of merger events. Addressing these various 
issues to improve the accuracy of the theoretical predic- 
tions is beyond the scope of this paper, but is certainly 
of interest for future work. 
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